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摘要 : 运动 目标 成 像 在 实际 应 用 中 具有 重要 作用 ， 而 如 何 获取 高 质量 运动 目标 图 像 是 该 领 


域 研究 中 的 一 个 热点 问题 。 本 文采 用 行 扫描 的 采样 方式 ， 通 过 构造 运动 测量 和 矩阵， 建立 一 


种 基于 压缩 感知 理论 的 运动 物体 成 像 模型 ， 并 通过 仿真 及 实验 ,验证 了 该 模型 对 于 恢复 运 


动物 体 图 像 信息 的 可 行 性 。 实 验 结果 证 明 ， 该 方法 可 获得 高 质量 的 运动 物体 成 像 。 通 过 引 


入 图 像 质 量 评价 标准 ,分 析 了 运动 物体 成 像 质 量 与 速度 之 间 的 关系 。 本 文 还 将 该 方法 与 普 


= 


通 压缩 感知 算法 进行 比较 ， 结 果 证 明 ， 在 相同 速度 下 ， 该 方法 的 成 像 质量 更 高 。 因 而 ， 该 


成 像 方法 在 无 人 机 对 地 观测 、 产 品 线 视频 监测 等 领域 有 着 很 好 的 应 用 前 景 。 


关键 词 : 压缩 感知 ， 运 动 目 标 成 像 ， 行 扫描 测量 矩阵 


1 引言 


压缩 感知 (Compressed Sensing, CS) 理论 是 由 Candss、Donoho、Tao 等 人 提出 的 一 


种 新 型 信号 重建 理论 。 该 理论 指出 ， 只 要 信号 是 可 压缩 的 或 在 某 个 变换 域 下 是 稀 玻 的 ， 


可 通过 随机 的 观测 矩阵 将 高 维 信号 投影 到 低 维 空间 上 , 并 通过 求解 一 个 优化 问题 从 少量 的 


投影 值 重 建 原始 信号 。 由 于 其 打破 了 夺 奎 斯 特 采 样 定理 的 限制 ,所 以 该 信号 恢复 理论 成 为 


近年 来 一 个 热点 科学 问题 5 1。 基于 CS 理论 的 典型 应 用 当 属 Rice 大 学 提出 的 单 像素 相机 模 


型 ， 即 仅 用 一 个 点 探测 器 收集 信息 便 可 成 像 和 "。 还 有 学 者 将 其 应 用 于 量子 成 像 ， 遥 感 成 


像 ， 医 学 影音 等 领域 …。 与 其 他 成 像 方式 一 样 ， 获 得 高 质量 成 像 也 是 CS 算法 追求 的 主要 


目标 之 一 。 然 而 ， 上 述 CS 成 像 研究 均 是 针对 静止 目标 而 言 。 


在 实际 成 像 中 ， 成像 目 标 大 多 数 是 运动 的 ， 因 而 获取 运动 目标 的 高 质量 成 像 具 有 重要 


意义 。 比 如 在 无 人 机 对 地 观测 、 产 品 线 视频 监测 等 过 程 中 ,成 像 系统 与 探测 目标 之 间 的 相 


对 运动 会 导致 获取 的 目标 成 像 质量 退化 ,产生 运动 模糊 ， 所 以 运动 目标 的 成 像 质 量 在 实际 


DOA PAA gg s f jg Ute 
近年 来 ， 不 少 学 者 提出 过 有 关于 CS ISR RRA, ECE tia Be E AE Se LH A 
用 散 斑 大 小 和 速度 检索 可 实现 运动 鬼 成 像 的 去 模糊 " 1, BIAS AC IIR or E EHE LB T 


光学 遥感 运动 补偿 压缩 成 像 ， 通 过 对 运动 条 件 下 压缩 采样 的 精确 建 模 ， 实 现 了 运动 补偿 压 


缩 成 像 ， 有 效 提高 了 光学 遥感 的 成 像 质量 "以 及 本 实验 室 研 究 小 组 提出 的 基于 互补 采 


样 的 运动 目标 压缩 成 像 "1。 传 统 CS 成 像 算法 普遍 采用 区 域 扫描 成 像 ， 但 是 如 果 待 重建 图 


像 像素 规模 较 大 时 ， 测 量 和 矩阵 规模 也 会 很 大 ， 则 在 存储 测量 矩阵 时 会 占用 系统 很 大 内 存 。 


本 文 提出 采用 采用 行 扫描 的 采样 方式 [23, 24], ， 通 过 构造 运动 测量 和 矩阵， 建立 一 种 新 的 基 


于 压缩 感知 的 运动 物体 成 像 模型 , 采用 行 扫描 采样 方式 不 仅 可 以 使 测量 矩阵 大 幅 减 少 存储 


量 并 提高 测量 采样 时 间 。 我 们 还 在 数值 模拟 的 基础 上 ， 引 入 数字 微 镜 器 件 记 ”7， 搭 建 了 一 
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套 原理 验证 实验 装置 ， 验 证 了 该 模型 可 用 于 重建 高 质量 的 运动 物体 图 像 。 
2 理论 


2.1 压缩 感知 理论 


压缩 感知 理论 (Compressed Sensing, CS) 作为 一 种 信息 处 理 技术 自 2006 年 被 提出 以 


来 就 被 有 效 地 应 用 于 许多 光学 领域 。 该 理论 指出 ， 对 稀 玻 或 可 压缩 的 信号 进行 少量 非 适应 


性 的 线性 测量 ， 所 获取 的 信息 就 足够 用 来 重建 信号 ， 这 突破 了 传统 的 奈奈 斯 特 采样 定理 。 


在 压缩 感知 理论 中 ， 通 过 利用 信号 的 稀 朴 性 先 验 知识 和 优化 算法 ， 可 以 通过 m 次 测量 恢 


复 n 维 信号 m<) ,也 就 是 说 可 以 实现 亚 采 样 重建 物体 信息 。 在 压缩 感知 数学 模型 中 ， 


一 个 一 维 信号 可 以 在 一 个 稀 朴 基 到 下 稀 玻 表示 : 


x= (1) 


HP x EX mA, MELER DANA: 


y-Ax-e (2) 


其 中 yy 为 m HEX8, AZg— mXn5BEE, BUM, Km 可 以 小 于 n，e 为 噪声 ， 当 


矩阵 4 满足 RIP (restricted isometry property) 条件 时 站， 信号 x 可 通过 下 面 的 优化 目标 
函数 得 到 : 


vin! raw + efx o 


其 中 z 为 一 常数 因子 ， 


言 号 的 黎 下 表示 是 压缩 感知 可 亚 采 样 有 效 恢 


复 物体 信息 的 重要 前 提 和 关键 点 , 此 外 采样 矩阵 的 设计 构造 和 重建 算法 也 是 压缩 感知 的 核 
心 内 容 " 亿 4。 压缩 感知 理论 凭借 其 独特 的 优势 在 光学 领域 得 到 了 成 功 的 应 用 1。 


2.2 基于 改进 行 扫描 测量 矩阵 的 运动 目标 压缩 成 像 


上 面 介绍 的 算法 已 经 广泛 使 用 在 对 静态 物体 的 压缩 感知 成 像 中 , 对 于 运动 物体 的 成 像 


需要 重新 定义 。 本 文 在 行 扫描 的 基础 上 建立 了 基于 CS 的 运动 物体 成 像 模型 。 模 型 中 引入 


参数 y Cpixels/sampling) "”"， 来 表示 运动 物体 的 速度 ， 并 规定 物体 在 某 一 直线 方向 上 运 


动 ， 行 扫描 方向 与 物体 运动 方向 为 垂直 关系 ， 如 图 〈1) 所 示 。 本 文 做 出 这 样 的 扫描 规定 


是 符合 机 载 扫 搬 或 产品 线 视频 监测 等 应 用 的 实际 情况 的 。 因 为 物体 连续 运动 ， 而 探测 器 不 


能 保证 完全 采集 到 运动 物体 的 全 部 信息 ， 根 据 图 像 积 分 理论 站， 假设 物体 在 dt 时 间 内 运 


Dn 


动 距离 为 4 ， 如 图 1 所 示 ， 则 速度 可 以 表示 为 公式 〈4) 所 示 : 


di 


d (4) 
设 图 像 每 行 像素 数 为 w ， 才 长 度 内 的 像素 数 为 必 , M: 
ds =s,-dl (5) 
将 公式 (4) 代入 上 式 有 ; Psy (6) 
将 上 式 两 边 同时 对 s 求 积分 ， 便 得 到 物体 直线 运动 的 连续 函数 ， 如 公式 (7) HR: 
s= [sovat (7) 
总 长 


每 行 像 素数 : 50 


图 1. 物体 运动 示意 图 


本 文 规定 竺 成像 的 物体 是 一 幅 像素 规模 为 nxXn BUEME Xx. BARAT HATES, Br 


以 在 此 x 每 次 仅 有 一 行 宽度 的 像素 内 容 参 与 成 像 ,经 过 多 次 行 扫描 过 程 ,最 终 拼 成 一 幅 nxn 


的 图 像 。 因 为 物体 在 连续 运动 ， 实 际 参与 每 次 行 扫 成 像 的 并 不 是 一 行 像素 内 容 ， 而 是 多 于 


一 行 ， 多 行 像素 内 容 混 琶 恢 复出 重建 图 的 一 行 像素 内 容 ， 如 果 按 照 正常 CS 成 像 直接 求解 


会 使 获取 的 物体 图 像 严 重 退化 。 


所 以 要 想 通 过 CS 算法 恢复 出 高 质量 的 运动 物体 成 像 ， 就 需要 找 出 参与 成 像 过程 中 行 


与 行 之 间 的 关联 ， 这 也 是 本 模型 采用 行 扫描 采样 方式 的 原因 之 一 。 实 际 采样 过 程 中 ， 行 扫 


描 的 每 次 调制 过 程 都 是 物体 的 两 行 像素 内 容 在 参与 成 像 ， 这 里 我 们 称 之 为 “当前 行 x,” 


和 “下 一 行 x, ”， 而 且 在 当前 行 x, 完全 运动 出 探测 器 的 扫描 区 域 之 前 ，x 和 xv 所 包含 的 


像素 内 容 大 小 一 直 在 变化 ， 设 5 表示 行 扫描 每 次 采样 过 程 中 探测 器 收集 到 的 像素 信息 ， 代 


入 公式 CD 可 将 5 表示 如 下 ， 


S=(s, - f s,vodt)+ | savodt (8) 
0 0 


RE s, Als, 分别 表示 当前 行 x, 和 下 一 行 x, 所 包含 的 全 部 像素 信息 ， 公 式 (8) 中 等 号 右边 


第 一 部 分 w - | sd 表示 物体 连续 运动 中 ， 探 测 器 收集 到 的 当前 行 参与 成 像 的 像素 信 


息 ， 即 用 当前 行 的 全 部 信息 s, 减 去 运动 出 去 的 像素 信息 (积分 部 分 ) 得 到 ;第 二 部 分 


[ssvwdt 表示 收集 到 的 下 一 行 x, 参与 成 像 的 像素 信息 ， 即 《〈 积 分 部 分 ) 表示 下 一 行进 入 到 


探测 器 探测 范围 的 像素 信息 。 将 公式 〈8) 对 应 到 CS 成 像 算法 上 表示 为 : 


1 1 1 1 
y - Al, — | x,vydt) + | xgvodt+e=AL(1— f vydt)x, + (| vodt)x,]+e (9) 
0 0 0 0 


即 一 次 采样 期 间 ， 和 采集 到 的 信息 为 RI oe, ASM BAA AI, AY WG IEE 
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X52], x, 与 xj 像素 信息 所 占 的 权 值 比例 一 直 在 变动 。 本 文中 物体 匀速 运动 ， 即 vy 7g 


数 ， 由 公式 “9) 中 的 积分 关系 ， 可 以 很 容易 推导 出 如 下 公式 : 


y-A(p:x,*q:x,)*e- (A: p) x, *CA:q? xı +e (10) 


其 中 p ，4 分 别 表示 行 扫描 的 每 次 采样 中 当前 行 x, 和 下 一 行 x 被 探测 器 收集 到 的 像素 内 


容 所 占 的 权 值 ， 即 公式 (9) 中 积分 计算 推导 出 的 1- f vedt A f vdt x, 完全 运动 出 探测 


0 0 


器 的 扫描 区 域 时 ， 这 时 完成 了 一 次 行 扫描 ， 恢 复出 一 行 图 像 信 息 ， 需 要 说 明 的 是 这 里 恢复 


出 来 的 一 行 图 像 为 X1 与 X2 拼接 的 一 行 图 像 ， 用 和 矩阵 可 表示 为 [x,%]， 即 像素 范围 是 1X 


2n, 在 呈现 成 像 结果 时 ， 将 其 转换 成 2xn 维度 ， 即 直接 恢复 出 前 两 行 图 像 信息 。 下 一 次 行 


扫描 时 ， 原 图 像 的 第 三 行 变 成 了 此 次 行 扫描 的 “当前 行 ”， 并 以 此 类 推 完成 m2 次 隔行 扫 


描 ， 即 可 恢复 出 完整 图 像 信息 。 采 样 到 最 后 一 行 时 ， 此 时 视野 区 域外 也 可 参与 计算 ， 呈 现 


NEA 
2 
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9 过程 写成 矩阵 形式 为 : 


y =(A-P)x,+(A-Q)x, +e =[A-P,A-Q][x,,x,]+e (11) 


根据 CS 理论 ， 我 们 把 公式 (11) 中 的 [4. 己 4.0] 称 为 改进 运动 测量 矩阵 ， 用 妃 表 示 , B: 


B -[A- P, A-Q] [PO] [A. A] (12) 
D oe Pi Q . qi 

其 中 ， P=| : >, 3 ; Q2: 7 3 (13) 
Pm pos Pn Gin b Gin 


根据 公式 (13) WA, Spe sump yy sin S3 REIR RS HARE P.O 值 的 设置 决定 ; 


FESS A BOR REA, OP, O 矩阵 进行 调整 即 可 ;， 其 中 p, dn 分 别 对 应 每 次 行 


AKE m 次 采样 过 程 中 x, M x 像素 内 容 所 占 的 权 值 , 在 当前 行 x, 完全 运动 出 探测 器 的 扫 
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描 区 域 之 前 ，xz 和 xs 所 占 的 权 值 一 直 在 变 ， 即 x, 参与 成 像 的 像素 一 直 在 变 少 ， 权 值 在 变 


小 ， 相 反 x, 的 权 值 一 直 在 变 大 。 


站 De Ge dg 

B=[P,Q] [A,A]=| : 7. i i 0 i D[A,A] (14) 
Pn iun Pn Qin Dd Qin 
x =[x, 2.7 (15) 


(16) 


AR 14) 中 的 4 为 普通 CS 算法 中 的 测量 矩阵 ， 和 矩阵 维度 是 mxn, B 为 本 文 模型 中 


的 改进 运动 测量 矩阵 ， 维 度 为 mx2n， 只 要 B 满足 相关 约束 条 件 ， 根 据 CS 重建 算法 ， 便 


可 重建 出 运动 目标 在 不 同 速度 下 的 图 像 。y 为 每 次 行 扫描 的 测量 值 。x 表示 由 x 和 x 组 


成 的 一 个 2nx1 维 的 列 向 量 ， 即 实验 中 重建 出 的 图 像 信 息 。e 为 加 入 噪声 ， 这 里 取 高 斯 噪 


声 ， 噪 声 标 准 差 是 观测 值 标准 差 的 5 多 ， 表 示 为 公式 〈17) 。 


Std (e) = 0.05Std(y) (17) 


以 上 的 改进 CS 算法 ,充分 利用 了 行 扫描 采样 方式 在 处 理 行 与 行 之 间 像素 关联 问题 的 优势 ， 


有 效 地 解决 了 运动 物体 成 像 问题 。 


3 实验 与 数据 分 析 


为 验证 本 文 提 出 的 基于 行 扫描 的 改进 测量 矩阵 运动 物体 成 像 模型 ， 利 用 DMD 数字 


微 镜 器 件 ，Digital Micro-mirror Device) 1 可 对 入 射 光 进 行 空 间 光 调制 ， 我 们 搭建 了 一 


套 基于 CS 的 运动 物体 成 像 实验 装置 ， 如 图 2 所 示 ， 采 用 讽 素 灯 作 为 光源 ， 光 源 经 过 衰减 


片 和 滤 光 片 后 通过 会 聚 透镜 将 光 会 聚 到 DMD 上 一 定 像素 区 域内 。 实 验 中 所 采用 的 DMD 


的 微 镜 阵 列 规模 是 1024X768， 微 镜 尺 寸 为 13.68hmx13.68hum， 每 个 微 镜 可 以 单独 控制 旋 


转 十 12" 方 向 或 者 一 12" 方 向 ， 把 这 两 种 状态 记 为 1 和 0， 分 别 对 应 “ 开 ” 和 “ 关 ” 状 态 。 


“ 开 ” 状 态 时 ，DMD 将 光 反 射 至 后 续 的 收 光 系统 中 ; “ 关 ” 状 态 时 ， 其 反射 区 不 被 收集 。 


实验 中 为 简化 分 析 ， 让 待 成 像 物 体 为 一 幅 二 值 虚拟 物体 图 像 ””， 如 图 3 Ca) 所 示 为 


一 大 写字 母 A， 像 素 规模 为 32x32。 每 次 行 扫 的 图 像 区 域 为 1X32 的 像素 范围 ， 实 验 中 采 


取 “ 像 素 扩 充 ” 方 法 ,将 待 成 像 物体 的 1x32 像素 区 域 等 比例 扩充 24 倍 对 应 到 DMD 上 24 


X768 的 像素 区 域 上 ， 区 域外 都 处 于 “ 关 ” 状 态 。 实 验 中 我 们 设 定 物体 运动 速度 为 


5; 0041667 pixels/sampling( 即 每 次 行 扫描 进行 24 次 采样 ， 这 期 间 物 体 运 动 了 一 个 像素 


宽度 的 距离 ) "“ ， 对 应 采样 率 为 对 应 到 扩充 之 后 的 DMD 上 ， 相 当 于 每 次 采样 时 间 内 ， 


物体 运动 了 一 个 DMD 微 镜 尺 寸 宽 度 的 距离 。 这 时 : 


| 23 23 | | 1 1 | 

Pals A Bj Q= 2-4 (18) 
0 0 24 24 
L24 24 [24 — 24] 


则 公式 (14) 的 改进 运动 测量 矩阵 可 表示 为 公式 〈19) 如 下 : 


B -[P,.Q, | A. A] (19) 


上 式 中 关于 测量 和 矩阵 A. 的 选择 一 条 最 基本 的 原则 是 满足 RIP 准则 , 根据 以 往 学 者 的 研究 成 


果 ”可 知 服从 伯 努 利 分 布 的 -171 BENLE Be i AE RIP 准则 并 且 容 易 实现 而 被 广泛 使 用 ， 


它 恰 好 能 够 满足 DMD 这 种 二 位 空间 光 调 制 器 的 控制 需求 。 另 外 ， 随 机 二 进 制 矩阵 的 元 素数 


值 仅 需 一 位 即 可 表示 ， 可 以 有 效 的 节约 存储 空间 和 运算 时 间 。 因此 本 文中 的 A 设计 为 


24x32 的 -1/1 ty NL AEB - 


再 用 B 对 每 行 待 成 像 物体 的 像素 内 容 进行 调制 得 到 Y,， 将 24 AY, XE EF DR E 
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DMD 上 ， 依 次 对 待 成像 物 体 的 每 一 行 像素 宽度 的 图 像 进行 压缩 传 感 测量 与 重建 ， 每 次 行 


扫 测 量 数 为 24， 总 共 32 行 ， 即 进行 16 次 隔行 扫描 ， 由 公式 〈19) 的 改进 测量 矩阵 可 看 


出 每 次 行 扫描 采样 率 为 24/64=0. 375。 实 验 中 经 过 DMD 之 后 的 反射 光 被 收集 透镜 L4 收集 ， 


利用 光电 桶 探测 器 探测 总 光 强 ， 如 图 2 所 示 。 


光源 滤 光 片 L， ”衰减 片 L。 SRL, DMD EMH) 


收集 透镜 La 


N PEE Q> 桶 探测 器 Dt 
SS 


A 


WW! 


2 实验 装置 图 


最 后 ， 测 量 值 y 和 改进 运动 测量 矩阵 B 作为 重建 算法 软件 的 输入 ， 利 用 TVAL3 算法 


求解， 重建 出 待 测 运动 图 像 信息 ， 将 每 行 恢复 出 的 像素 信息 拼接 ， 便 得 到 完整 恢复 图 


像 ， 如 图 3 (c) 所 示 。 此 外 我 们 还 在 不 同 速 度 下 比较 了 普通 CS 算法 和 本 文 改进 CS 算法 


对 于 恢复 运动 物体 图 像 的 能 力 ， 实 验 表明 ， 在 不 同 速度 下 改进 CS 算法 均 具 有 更 好 的 运动 


物体 重建 能 力 ， 如 图 3 所 示 。 


图 3 实验 重建 图 


chinaXiv:201703.00268v1 


Ca) 原 图 ， 像 素 规模 32x32; (b) 普通 CS 算法 下 速度 为 1/24 pixels/sampling (采样 率 
为 0.3750) 的 重建 图 ， (c) 改进 算法 下 速度 为 1/24 pixels/sampling (采样 率 为 0.3750) 
的 重建 图 ， (Cd) 普通 CS 算法 下 速度 为 /12 pixels/sampling (采样 率 为 0.0833) 的 重建 
Al; (e) 改进 算法 下 速度 为 /12 pixels/sampling (采样 率 为 0.0833) 的 重建 图 . 


实验 证 明 ， 根 据 本 文 提供 的 基于 行 扫描 的 改进 测量 矩阵 CS 成 像 方法 ， 在 一 定 速度 内 


可 以 实现 高 质量 的 运动 物体 图 像 。 此 外 我 们 还 通过 仿真 实验 验证 了 该 改进 算法 的 有 效 性 。 
仿真 实验 中 我 们 选取 了 经 典 Lene 灰 度 测试 图 其 中 一 部 分 的 120X 120 范围 作为 成 像 目 标 ， 
并 设 定 一 组 匀速 运动 的 速度 值 ， 对 普通 CS 算法 和 改进 算法 进行 比较 ， 得 到 如 图 4 所 示 的 


运动 物体 重建 图 。 可 见 本 文 提出 的 改进 算法 对 于 重建 运动 目标 成 像 具 有 显著 效果 。 


(b1) (b2) (cl) (c2) (d1) (d2) 


图 4. 仿 真 重建 图 


(a) 原始 图 像 64X64; (b1) 改进 CS 算法 下 速度 0.01163pixels/sampling 重建 图 ;(b2) 


普通 CS 算法 下 速度 0.01163pixels/sampling 重建 图 ， (cl) 改进 CS 算法 下 速度 
0.01724pixels/sampling 重建 图 ;，(c2) 普通 算法 下 速度 0. 01724pixels/sampling 重建 
Al; (d1) 改进 CS 算法 下 速度 0. 02632pixels/sampling 重建 图 ;，(d2) 普通 算法 下 速度 


0. 02632pixels/sampling 重建 图 . 
4 性 能 对 比 
为 了 有 效 评价 重建 图 像 质量 的 好 坏 , 本 文 还 引入 了 均 方差 (Mean Square Error,MSE) , 


定义 式 为 公式 〈20) ， 这 是 基于 像素 点 间 误 差 敏感 的 全 参考 图 像 质量 评价 指标 ， 即 把 像素 


误差 看 作 是 图 像 质量 的 亏损 。MSE 值 越 小 ， 重 建 的 图 像 质量 越 好 。 广 ， 


1 M N i 
MSE — v 222,008 tia) (20) 


=l k=l 


AP: x 表示 原始 图 像 信号 ，x 表示 重建 图 信号 ， 图 像 大 小 为 MxN， jo 分 别 代表 图 像 


的 行 和 列 ， 分 析 结 果 如 图 5 所 示 。 给 出 普通 CS 算法 和 本 文 改进 运动 测量 年 阵 算法 的 MSE 
值 随 速 度 变化 曲线 ， 随 着 速度 提高 ， 图 像 重 建 质量 均 有 所 减弱 ， 而 本 文 改 进 运动 测量 矩阵 


算法 的 MSE 值 明 显 一 直 低 于 普通 CS 算法 。 所 以 ， 本 文 改进 算法 在 提高 运动 物体 CS 成 像 


质量 上 有 显著 成 效 。 同 时 从 图 5 中 还 可 以 看 出 当 速 度 逐 渐变 大 时 , 图 像 重建 质量 越 来 越 差 ， 


这 是 因为 当 速度 逐渐 变 大 时 ， 每 次 行 扫描 所 采样 到 的 每 一 行 的 像素 内 容 越 来 越 少 ， 导 致 重 


建 图 像 的 退化 越 来 越 严 重 。 从 图 G) 可 看 出 当 运 动 速度 低 于 一 定 值 时 ，MSE 值 随 速度 大 


小 的 变化 缓慢 ， 这 是 因为 物体 运动 速度 慢 时 相当 于 行 扫描 采样 数 已 经 足够 CS 算法 的 重建 


需求 ， 同 样 当 速度 快 于 一 定 值 时 ， 可 看 出 ME 曲线 变化 缓慢 ， 这 是 因为 当 速 度 快 于 一 定 值 


以 后 ， 每 次 行 扫描 的 采样 数 无 限 少 ， 此 时 的 重建 图 已 经 接近 一 幅 随 机 图 ， 所 以 速度 再 快 对 


于 MSE 值 的 影响 也 很 微弱 了 。 可 见 随 着 速度 的 加 快 ， 采样 数 势 必 越 来 越 少 ， 采 样 率 越 来 越 


低 ， 重 建 质量 也 逐渐 变 差 。 


一 Se 一 改进 算法 


MSE 


0 | | | | | 
0.01 0.012 0.014 0.016 0.018 0.02 0.022 0.024 0.026 0.028 


速度 〈pixels/sampling ) 
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图 5. MSE 随 速 度 变化 曲线 
5 FARER RR 


以 上 描述 均 已 匀速 直线 运动 物体 为 前 提 , 本 节 将 本 文 提 出 的 模型 用 于 非 匀 速 运动 直线 


运动 成 像 的 计算 ， 并 给 出 仿真 实验 结果 。 由 公式 〈8) 可 以 得 出 ， 非 匀速 运动 时 将 公式 中 


的 常量 vO 写成 随时 间 变 化 的 变量 vt 即 可 ， 如 公式 (21) 所 示 。 


S=(s, - [5,40 + svar (21) 
0 0 


2275 EX ^D Bde s MA OR Th AE AY UA HE Se SIE. 只 要 各 个 时 


刻 的 速度 己 知 ， 对 公式 〈13) 中 P, Q 矩阵 进行 调整 即 可 得 到 采样 过 程 中 描述 当前 行 与 下 一 


行 像素 内 容 所 占 权 值 关 系 的 矩阵 。 本文 给 出 加 速度 已 知 的 匀 变 速 直 线 运 动物 体 成 像 的 计算 


过 程 和 仿真 结果 。 设 加 速度 为 a (单位 为 bixels/sampling ) ， 初 速度 为 0， 则 任 一 时 刻 


的 瞬时 速度 可 表示 为 公式 (22) 所 示 : ”v=at (22) 


这 里 的 t 应 理解 为 采样 的 次 数 。 


将 公式 (22) WAAR (21) 并 求解 积分 可 以 得 到 S — s 0-7 a) es, rar (23) 


同样 s,s, 还 是 分 别 代表 当前 行 和 下 一 行 所 包含 的 全 部 像素 信息 , 则 同 匀 速 运动 一 样 ， 


对 应 到 CS 成 像 算 法 上 表示 为 


y-A(p,:x,*q,:x,)*-e-(A: p) x,-(A:q) x, +e=[A-P,A-O][x,,x,]+e (24) 


根据 上 式 可 以 看 出 ， 非 勾 速 运动 时 的 计算 过 程 和 义 速 运动 成 像 的 计算 过 程 基 本 相似 ， 


加 速度 已 知 ， 则 任 一 时 刻 的 速度 也 可 得 到 ， 只 要 对 当前 行 和 下 一 行 所 占 像 素 的 权 值 进行 相 


应 调整 即 可 ， 也 就 是 对 矩阵 Pt, Qt 中 的 元 素 乘 以 相应 比例 系数 进行 调整 。 本 文 给 出 一 系列 


加 速度 值 下 的 匀 变 速 直线 运动 物体 成 像 的 仿真 结果 。 
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图 6. 5J4EXE ELA oh Hj EA 


Ca) JAI, 160X160; (b) a=0.027027 重建 图 ; (c) a=0.030303 重建 图 ;(d) a=0.034483 重建 图 ; 


(e) a=0.076923 重建 图 . 


6 结论 与 展望 


本 文采 用 行 扫 描 的 采样 方式 ,通过 构造 运动 测量 矩阵 ,建立 了 基于 压缩 感知 的 运动 物 


体 成 像 模 型 。 并 搭建 了 基于 DMD 的 原理 验证 实验 光路 系统 ， 验 证 了 该 模型 在 重建 运动 物 


体 成 像 方 面相 较 普通 CS 算法 具有 更 高 恢复 能 力 。 同 时 也 发 现 当 速度 快 于 一 定 值 时 ， 重 建 


图 像 质量 发 生 较 显 著 的 退化 现象 ， 而 且 此 模型 目前 分 别 对 匀速 水 平 、 垂 直 运 动 和 勾 变 速水 


平 、 重 直方 向 运动 物体 成 像 的 重建 有 效 ， 若 物体 治 倾斜 方向 运动 ， 需 要 对 模型 进行 相应 修 
正 ， 针 对 以 上 缺陷 ,我 们 将 在 今后 的 工作 中 展开 进一步 研究 。 鉴 于 本 文 提 出 的 基于 行 扫描 


的 改进 测量 和 矩阵 压缩 感知 运动 成 像 算 法 的 优势 , 有 理由 相信 将 在 实际 成 像 领域 中 得 到 广泛 
MHo 
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Moving target compressive imaging based on 


improved row scanning measurement matrices 
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(1. Key Laboratory of Electronics and Information Technology for Space Systems, 
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2. University of Chinese Academy of Sciences, Beijing 100190, China; 


3. School of Physics, Beijing Institute of Technology, Beijing 100081, China) 


Abstract: Moving target imaging (MTI) plays an important role in practical applications, and 


how to capture high quality dynamic images of the target becomes a front-burner issue in this 
field. In this paper, a new MTI model based on compressed sensing is proposed to improve its 
reconstruction quality, applying a sampling scheme of the row-scanning together with a motion 
measurement matrix constructed by our own. The simulation and the experiments show that a 
relatively higher quality of the moving target imaging can be achieved through this approach. 
Furthermore, an evaluation criterion of reconstructed images is introduced to analyze the 
relationship between the imaging quality and the moving speed of the target. By comparison, it 
can be concluded that the performance of our algorithm is better than that of traditional 
compressed sensing algorithms under the same moving speed condition. Therefore, the results 
suggest that our imaging method has a great application prospect in the earth observation of 


unmanned aerial vehicles, video monitoring in the product line and other fields. 


Keywords: compressed sensing; moving target imaging; row-scanning; 


measurement matrix 
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